Slow dynamics in cylindrically confined colloidal suspensions 
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Abstract. We study bidisperse colloidal suspensions confined within glass microcapillary tubes to model the glass transition 
in confined cylindrical geometries. We use high speed three-dimensional confocal microscopy to observe particle motions 
for a wide range of volume fractions and tube radii. Holding volume fraction constant, we find that particles move slower in 
thinner tubes. The tube walls induce a gradient in particle mobility: particles move substantially slower near the walls. This 
suggests that the confinement-induced glassiness may be due to an interfacial effect. 
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1. INTRODUCTION 

Understanding the glass transition is one of the endur- 
ing questions of solid-state physics [1-6]. The problem is 
simply stated: in some cases, when a hot viscous liquid 
is cooled, the viscosity rises dramatically but smoothly 
as a function of temperature. At some temperature the 
viscosity is so large that the sample appears like a solid; 
this identifies the glass transition temperature. The same 
phenomenon can likewise be induced by increasing the 
density (increasing the pressure) [7]. In contrast to a reg- 
ular phase transition which occurs at well-defined tem- 
peratures and pressures, the glass transition can depend 
on details such as the cooling rate. Likewise, while phase 
transitions are signaled by abrupt changes in the sample 
properties or their derivatives, the properties of a glass- 
forming material such as viscosity and diffusivity change 
smoothly, and to an extent the definition of the transition 
temperature (or pressure) is a bit arbitrary [8, 9]. One 
of the key questions is what is changing microscopically 
that is responsible for the macroscopic changes in vis- 
cosity; no structural length scale has yet been found that 
would clearly explain the viscosity change [10, 11]. 

One clever way to probe length scales is to confine 
a sample: rather than studying a macroscopically large 
sample (a "bulk" sample), studying a microscopic-scale 
sample. Many experiments show that glasses change 
their properties when their size is sufficiently small, both 
for small-molecule glasses and polymer glasses [8, 12- 
14]. One of the key observations is that the glass tran- 
sition temperature T g changes for confined samples. In 
some cases T g increases: confined samples are glassier. 
However, in other cases T g decreases. The key differ- 
ence explaining the increase or decrease seems to be the 
boundary conditions [14]. Samples with free surfaces, 
such as thin free-standing polymer films, are less glassy 
(lower T g ). Samples confined to pores or on substrates 
can be more glassy (larger T g ), in particular if the sample 



molecules form strong chemical bonds to the confining 
boundaries. 

We wish to use colloidal samples as a glass-forming 
system which can be studied in confinement. Colloidal 
suspensions are composed of solid particles in a liquid. 
As the particle concentration is increased, the sample be- 
comes more and more viscous [15-18]. Above a critical 
concentration, the sample behaves as a glass, and a large 
number of similarities have been observed between the 
colloidal glass transition and glass transitions of poly- 
mers and small molecules [9]. The most widely studied 
colloidal glass transition is that of hard-sphere-like col- 
loids, and the control parameter is the volume fraction 
[19]. The glass transition point has been identified as 
fy g ~ 0.58, with simulations demonstrating that this re- 
quires some polydispersity [20, 21]. 

Experimentally, the colloidal glass transition shifts to 
lower volume fractions in confined samples: confinement 
makes colloidal samples glassier [22-27]. This has been 
studied exclusively in parallel-plate geometries, where 
samples are confined between two glass walls that are 
closely spaced. Often these experiments use bidisperse 
samples (mixtures of two particle sizes) so that the flat 
walls do not induce crystallization [22, 24]. An alternate 
approach is to roughen the walls [25]. 

The geometry of these prior colloidal experiments 
most closely resembles thin films, which are used to 
study the glass transition of polymers or small molecule 
glass formers in thin slits. However, small molecule 
glass formers are more commonly studied by using 
nanoporous substrates; a variety of these nanoporous 
substrates are reviewed in Ref. [8]. Some of these sub- 
strates are quite disordered with pores of a variety of 
shapes and sizes. Others are ordered: for example, porous 
oxide ceramics have a regular lattice of cylindrical 
nanopores with well defined sizes [8, 28], as do anodized 
aluminum oxide membranes [29]. Experiments find that 
confinement in cylindrical pores can both enhance or di- 



minish glassy behavior [28]. Simulations show that the 
boundaries play an important role in this: rough walls 
that frustrate layering of particles result in glassier dy- 
namics, and smooth walls result in less glassy dynamics 
[30,31]. 

In this paper, we present a study of colloidal sam- 
ples confined in cylindrical glass tubes, to mimic the 
geometry of cylindrical nanopores. We use confocal mi- 
croscopy to observe both the structure and dynamics of 
the samples. Similar to prior colloidal work, we find 
that confined colloidal samples are glassier. In particu- 
lar, particle motion is dramatically slower at the capillary 
tube walls, demonstrating that we see an interfacial ef- 
fect. We use a bidisperse sample to prevent confinement- 
induced crystallization or other ordering, known to occur 
for monodisperse cylindrically confined spheres [32-37] . 
Nonetheless, the particles layer against the walls, and 
these layers slightly influence the motion in ways sim- 
ilar to previous observations [22, 24]. The data add to 
the analogy with confined small-molecule glass-formers. 
Additionally, they are of interest for colloidal suspen- 
sions themselves: the implication is that for microfluidic 
applications, it will be more difficult than anticipated to 
flow dense colloidal suspensions, as they will be glassier 
in small tubes than an equivalent bulk sample. Colloids 
confined in cylinders or small channels have been stud- 
ied before, but only in dilute concentrations [38] and/or 
in extremely thin channels that are only one particle di- 
ameter across [39, 40]. 



2. EXPERIMENTAL METHODS 

Our goal is to use hard-sphere-like colloids as a model 
system. Colloids have proved to be effective models with 
similarities to hard-sphere computer simulations; see [9] 
for a discussion. A hard-sphere system means that the 
particles do not interact with one another beyond their 
radius and are infinitely repulsive at contact [41-43]. 
An advantage is that the particle size can be selected to 
be ~ 1 jiim in radius: small enough to undergo random 
Brownian motion, yet still large enough to be imaged 
using microscopy [9]. 

We use poly(methyl methacrylate) (PMMA) spheres 
coated with a polymer brush layer that sterically sta- 
bilizes the particles, preventing them from aggregating 
[44]. We use a bidisperse mixture with particles of two 
different radii, large particles with radius a L = 1 .08 jxm 
and small particles with radius a s = 0.532 |im, helping 
us avoid crystallization [37]. The particles have a poly- 
dispersity of approximately 6% and additionally their 
mean radii and as are each uncertain by 1%. Our par- 
ticles are fluorescently dyed so that we can observe their 
motion with confocal microscopy. It is desirable to re- 
duce the influence of gravity in a colloidal suspension by 



density matching the colloid with the surrounding fluid. 
To do this we use a standard mixture of 85% (weight) 
cyclohexyl bromide and 15% decahydronaphthalene (de- 
calin, mixture of cis- and trans-). This solvent mixture 
also matches the particles' refractive index, which is nec- 
essary for the microscopy. To reduce the influence of 
electrostatic repulsion between the particles, we satu- 
rate the solvent mixture with tetrabutylammonium bro- 
mide (Aldrich, 98%) with a resulting concentration of 
~190juM[45]. 





FIGURE 1. Top: Composite photograph of capillary tube. 
The scale bar is 500 flm. Bottom: 2D confocal image from 
experiment 55, with volume fraction tot = 0.49. The scale bar 
is 10 /xm. 

Our sample chambers are glass capillary tubes as 
shown in Fig. l(top). The capillary tip was made using 
an automated pipette puller. To help the capillary tube fit 
onto a microscope slide, the large end of the capillary tip 
was cut off while ensuring not to break the thin end of the 
capillary tube. The thin end was often too thin, so it too 
was cut, leaving an opening a few microns in diameter. 
The shortened capillary tube was then dipped into a vial 
containing the colloidal suspension for 10 seconds or so, 
and the sample flows into the glass tube due to capillary 
forces. We use quick drying UV epoxy (Norland 81) to 
seal both ends of the tube and also to glue the tube to 
the slide. In fact, it is useful to cover the tip entirely with 
glue to ensure stability. The microscopy is not affected 
too much, as the sample and epoxy have similar indices 
of refraction ( 1 .495 and 1 .56 respectively). The result is a 



confined sample, such as shown in Fig 1 (bottom), which 
can be studied at different locations to give different tube 
sizes. The capillary tubes vary slightly in radius as a 
function of length (slope of about 1°), and the slight ta- 
per does not seem to affect our overall results. The taper 
is slight enough that we cannot observe it in any of our 
confocal images. Due to the pipette puller protocol we 
used, the capillary tubes sag slightly, and so their cross- 
section is slightly elliptical rather than circular. This does 
not seem to influence our results and will be discussed 
further below. When the tubes are filled with the sam- 
ple, some of the particles stick irreversibly to the tube 
walls due to van der Waals forces. This typically resulted 
in one complete layer of particles (of both sizes) coating 
the wall. During the course of our experiments, we do 
not observe any new particle sticking to the walls, nor do 
we observe any stuck particles becoming unstuck. 

We use a confocal microscope (Visitech vt-Eye) to 
study our samples. In a confocal microscope, laser light 
is scanned across the fluorescent sample and excites the 
dye to emit a different color light. The emitted light 
passes through a pinhole to remove the out-of-focus light 
and then is measured by a detector, a photomultiplier 
tube in our microscope. The sample is quickly scanned 
in x and y to acquire a two-dimensional (2D) image. 
Then the microscope focus is adjusted to scan more 2D 
images at different depths z, thus building up a 3D im- 
age. Being able to create 3D images makes confocal mi- 
croscopy a powerful tool for studying particle dynamics 
in the system. Each 3D scan takes about 2-3 seconds, and 
we typically take movies comprised of 400 of these 3D 
images. The particles are tracked in 3D using standard 
tracking techniques [46, 47]. Within each confocal im- 
age, the small and large particles are easy to tell apart 
[see Fig. 1 (bottom)] and so we can distinguish between 
them in our data [48]. 

We have two control parameters, the tube radius and 
the volume fraction. Complete descriptions of how these 
are measured are given in the Appendix; here, we sum- 
marize the key points. 

Our first control parameter is the tube radius. The po- 
sitions of the particles give us an accurate idea where 
the tube surface is. However, the tubes have an ellipti- 
cal cross-section rather than a circular cross section. We 
measure the major and minor axes R max and R m i n for 
each experiment. The ratio /? m ax/^min ranges from 1.14 
to 1 .39, with mean 1 .24. Because we are concerned with 
confinement effects, we report our data in terms of R m [ n 
in general, although both radii are listed for all exper- 
iments in Table I. Note that we report the radii corre- 
sponding to the maximum positions of the observed par- 
ticle centers: in general the particles at these maximum 
positions are the small ones (whose centers can get closer 
to the tube walls) and so the true tube sizes are larger by 
a s = 0.532 jum. 



TABLE 1. List of experiments, ordered by total vol- 
ume fraction 0tot- The volume fraction of the small 
particles can be determined by <j> s = 0t ot [l + ((1 — 

/)/ 7)W«s) 3 r'. where / = NsmallMot- The volume 
fraction of the large particles is then (j>i = to t — <Ps- The 
tube sizes R m i n and R mdx correspond to the maximum 
radii that the particle centers can reach; the physical 
tube walls are a distance « as further away. The value 
of (Ax 2 ) given is for At = 100 s, and corresponds to the 
information plotted in Figs. 2, 3. 



Expt 0tot fimin #max (Ax 2 ) N smaa /N M 
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Our other key control parameter is the volume frac- 
tion <j). We measure this in each data set by counting the 
numbers of small and large particles observed within a 
subvolume of the tube, and converting this to the volume 
fraction using the known particle sizes. Note that 1% un- 
certainties in the particle radii translate to 3% uncertain- 
ties of the volume fraction, and since each particle radius 
is uncertain, we have an overall systematic volume frac- 
tion uncertainty of at least 5% [49]. 

In summary, our experimental method is to image 
different portions of the same tube in hopes to get a 
constant volume fraction with differing tube radii, and 
to study different tubes with different volume fractions 
to understand the role of volume fraction. In practice, 
we determine these variables when the data are post- 
processed, and report the measured values in Table I. 

3. RESULTS 
3.1. Motion slows in confined samples 

By following the motion of all of the particles in 3D, 
we can observe how the motion depends on confine- 
ment. We quantify this by calculating the mean square 
displacement, defined as 

(Az 2 ) = ([z( f +A f )-z( f )] 2 ) (1) 

where (Az 2 ) is a function of the lag time At, and the 
angle brackets indicate an average over all times t and 



all particles. Here we are considering the z direction to 
be along the axis of the tube, primarily because this axis 
is perpendicular to the optical axis of the microscope 
and therefore has less position uncertainty 1 . The data are 
plotted in Fig. 2, where each family of curves correspond 
to a different volume fraction. Within each family, the 
slower curves correspond to narrower tubes: confinement 
induces glassier behavior. To further quantify this, we 
consider the specific value of (Ax 2 ) at a lag time At — 
100 s (which is arbitrary but chosen to match prior work 
[22]). This value is plotted as a function of minimum tube 
radius R m i„ in Fig. 3. The different symbols correspond 
to different ranges of volume fractions, and in general 
within each range, smaller tubes correspond to slower 
motion. Plotting the data against R max , (/? m ; n +/?max)/2, 
or V^min^max does not change the overall appearance of 
this graph significantly. It is intriguing to note that the 
magnitude of the effect is less significant than was seen 
in similar experiments with parallel plates [22]. Here, the 
strongest influence of confinement is to lower mobility 
by a factor of <~ 3 (for circles, squares, and triangles in 
Fig. 3). With parallel plates, mobility was lowered by a 
factor of ~ 40 for data with <j> w 0.46 [22], a volume 
fraction in the same range as our cylindrical data. We are 
unsure why the two experiments differ in the magnitude 
of mobility reduction. 
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FIGURE 2. Mean square displacements in the z direction 
(along the cylinder axis). Each family of curves is labeled at 
left by the volume fraction, and at right by the values of the 
minimum cylinder radius (in flm). For each curve, the volume 
fraction is within 0.02 of the labeled value. 



1 The mean square displacement for the other components are qualita- 
tively similar, except that they have more uncertainty which artificially 
increases the data at small lag times; see [49] for a discussion. 









8° ° : 










1.0 






□ <P o O -: 




: □ 














0.1 


: i . , 


A 
i , , 


i , , , i , , , i , , , i , , , i ; 




6 


8 


10 12 14 16 18 



Kin C« m ) 

FIGURE 3. The value of (Ax 2 ) at the time scale At = 100 s, 
plotted as a function of the cylinder minimum radius R m i n . 
The symbols indicate the volume fraction: circles are (j) = 
0.21 ± 0.02, squares are = 0.45 ± 0.02, diamonds are <j> = 
0.49 ±0.01, and triangles are = 0.53 ± 0.02. The circles 
correspond to the dotted lines in Fig. 2, the squares correspond 
to the dashed lines, and the triangles correspond to the solid 
lines. See Table I for a full listing of all volume fractions, 
mobility values, and other details. 

3.2. Motion is slower near walls 

Microscopy allows us to spatially resolve details of the 
motion. If confinement-induced slowing is a finite size 
effect, then the motion might be spatially homogeneous: 
the whole sample feels that it is small [14, 50]. If instead 
the confinement-induced slowing is an interfacial effect 
(due to the sample-wall interface), then particle motion 
would depend on where each particle is relative to the 
boundary. Of course, both effects could be present simul- 
taneously [28]. We check this by plotting the particle mo- 
bility (Ar 2 ) as a function of the distance s to the nearest 
wall in Fig. 4(a). It is immediately apparent that particles 
move slower when they are close to the wall (s — > 0), and 
a plateau value for the mobility isn't reached until sev- 
eral particle radii into the sample. Note that we are us- 
ing the full 3D mobility (Ar 2 ) = (Ax 2 + Ay 2 + Az 2 ), for a 
fixed lag time At = 30 s (for which we have more statis- 
tics), and now the angle brackets indicate an average over 
those particles with the specific value of s. The value of 
s is based on the initial position of the particle, at time t 
rather than t +At. It is probable that the mobility near the 
wall is very slightly enhanced, due to particles which dif- 
fuse away from the wall during At and thus enhance their 
mobility [38] . The different curves are for small and large 
particles as indicated. Not surprisingly, the smaller parti- 
cles are more mobile than the large ones. Our directly ob- 
served gradient in mobility is similar to that inferred from 
experiments on small molecule glasses [50-52], polymer 
glasses [53], and seen directly in simulations [30, 31, 54]. 
Our data show that the mobility changes over a distance 



of several particle diameters. 
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tangential to the tube wall (0 direction), and in the ra- 
dial direction (s direction). These components of mobil- 
ity are shown in Fig. 4(b). Here it is apparent that the ra- 
dial mobility is the most influenced by the oscillations of 

^large 

(s). Again, this is the same result seen before with 
fiat walls [22]. The particles at the maxima of ni arge (i) 
appear to be at favorable positions and their mobility is 
reduced, whereas those at the minima are in less favor- 
able positions with higher mobility. Motions in the z and 
9 directions are along contours of constant mean n and 
do not fluctuate with n{s). 

All of the results of Fig. 4 are qualitatively repli- 
cated in Fig. 5, which is data from a smaller radius 
tube (R m \ B = 7.8 jJ.m compared to R m i„ = 17.1 fim) and 
volume fraction only slightly larger than that of Fig. 4 
(0tot = 0.54 compared to 0.50). In Fig. 5(a), the mobility 
is lower near the boundary, lower for large particles, and 
oscillates with higher mobility corresponding to minima 
of niarge(s)- The data for the components [Fig. 5(b)] are 
noisier due to less statistics in the smaller tube, but again 
(Ar 2 ) shows a stronger anticorrelation with «i ar ge(' s )- The 
oscillations of «iarge(s) in Fig. 5(c) are more complex 
than those seen in Fig. 4(c), probably due to packing con- 
straints of a smaller tube. 



4. CONCLUSIONS 



FIGURE 4. (Color online) (a) Mobility as a function of dis- 
tance s from the wall, for the small particles and large particles 
as indicated, (b) Mobility for the large particles only, for the 
components of motion as indicated, (c) Number densities of 
the small and large particles. Note that the number density of 
the small particles has been multiplied by four. In all panels, the 
vertical dashed lines correspond to the local minima of nwe- 
They are spaced approximately 2.04 /im apart. The data are 
from experiment 51 with <p tot = 0.50 and R m [ n = 17.1 flm; see 
Table I for other details. 

A second feature of the data of Fig. 4(a) is that the 
curves oscillate. The oscillations are related to the fluc- 
tuations of the particle number density, as seen by com- 
paring Fig. 4(a) to Fig. 4(c). The latter shows the fluc- 
tuations of the number density of large and small par- 
ticles (solid and dotted lines, respectively). The mobility 
data in Fig. 4(a) are antic orrelated with niarge ( s ) '■ the local 
minima of n\ arge (s) are indicated by the vertical dashed 
lines, and correspond to local maxima of (Ar 2 ). The in- 
fluence of n sma ii is n °t seen, probably because the num- 
ber fraction of small particles is much less than that of 
the large particles (-/V sma ii /Marge = 0.07 for these data). 
The anticorrelation between number density and mobil- 
ity matches what has been seen in prior work on confined 
samples [22, 24]. 

Further insight into the mobility is found by splitting 
(Ar 2 ) into components parallel to the tube (z direction), 



Our results - in particular Figs. 4(a) and 5(a) - suggest 
that the slower motion of confined colloidal samples is 
due to an interfacial effect, where particles near the sam- 
ple walls are slowed. This agrees with a prior observa- 
tion that colloidal particle motion is slower near rougher 
walls [23], a result demonstrating that the nature of the 
confining walls plays a role and not merely the finite size 
of the sample chamber. It is unlikely this is merely a hy- 
drodynamic effect, as the magnitude of such an effect 
would only be a factor of ~ 2 — 4 in mobility and would 
not depend on R m i n [24, 38]. 

The clear gradient seen in Figs. 4(a) and 5(a) was not 
seen in prior work where colloids were confined between 
parallel walls [22]. It may be that the influence of con- 
finement is stronger in the cylindrically confined case: 
between parallel walls, there are two unconfined direc- 
tions, whereas in the cylinder there is only one uncon- 
fined direction [55]. Another possibility is that the cur- 
vature of the cylindrical walls introduces an effect not 
present with flat walls 2 . These geometrical differences 
are the only major differences between the experiments 



2 Although the hydrodynamic effect of the wall should be the same for 
flat or curved walls, sufficiently close to the walls. Prior simulations 
found essentially no difference in the diffusion constant between flat 
and curved walls, for particles within 1.5 diameters of the wall [38]. 
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FIGURE 5. (Color online) (a) Mobility as a function of dis- 
tance s from the wall, for the small particles and large particles 
as indicated, (b) Mobility for the large particles only, for the 
components of motion as indicated, (c) Number densities of 
the small and large particles. In all panels, the vertical dashed 
lines correspond to some of the local minima of n\ m „ e . The data 
are from experiment 64 with 0t ot = 0.54 and R m \ n = 7.8 )±m\ 
see Table I for other details. 



in this paper and those published earlier. As noted in 
Sec. 3.1, the other observed difference is that the slow- 
ing of motion in cylindrically confined samples (Fig. 3) 
is less pronounced compared to the observations between 
parallel plates [22]. This is counterintuitive given the ar- 
gument above that cylindrical confinement is "stronger" 
than parallel-plate confinement. This trend is the oppo- 
site of that seen for small molecule liquids [55]. Future 
experiments may be able to elaborate on the question for 
colloids: at one extreme, particle motion could be studied 
in a half-infinite system near a wall, whereas at the other 
extreme, particle motion could be studied in a spherical 
pore. The data presented in this paper suggest that chang- 
ing the dimensionality of the confinement in this way can 
result in interesting and qualitatively distinct behavior, 
in other words, enhancing a mobility gradient near walls 
while diminishing the overall confinement effect. 

Our results also imply that flowing colloidal suspen- 
sions through small cylindrical tubes will be harder for 



smaller tube radii. Additionally, the decrease in mobility 
near the walls perhaps implies an increase in the appar- 
ent viscosity of the sample near the walls, thus modifying 
the flow velocity profile in a nontrivial fashion. 

Note that we do not see any quantization effects: we 
do not see any particular change in the dynamics at 
any special ratios of particle sizes to tube sizes. This 
is in contrast to some theoretical predictions [56-58]. 
However, our data are only at a limited number of tube 
sizes, as shown in Fig. 3; our tubes are elliptical in cross 
section and so the ratio of particle size to tube size is not 
a constant for any given data set; and it is likely that such 
quantization effects are more subtle than we would be 
able to see in an experiment. 
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APPENDIX 

We explain in more detail how our experimental parame- 
ters are measured. In each case, we prepare samples and 
take data, measuring the tube radii and volume fraction 
from the data. 

Tube radius: As noted in Sec. 2, the tubes do not 
typically have a circular cross-section. Also, in general, 
the images of the tube are not precisely aligned with the 
xyz laboratory reference frame. To determine the radius, 
the position data are first rotated by — 4° around the 
x and y axes as necessary so that the z axis of the data 
corresponds with the tube axis. (Note that the z axis is 
not the optical axis of the microscope; rather, the z axis 
is within a few degrees of perpendicular to the optical 
axis.) Then the data are projected onto the xy plane 
and their center of mass is found. The x,y coordinates 
are converted to r, 9 and the maximum r is found as a 
function of 9. Following the procedure of Eral et al, 
in order to smooth the measured contour we fit r(9) to 
a Fourier series up to m = 4 modes [38]. The function 
r(9) is found to be well-described by an ellipse in all 
cases, and accordingly it is easy to determine the major 
and minor axes. 

Volume fraction: For our experiments, we take 
movies from different portions of several tubes. Ideally 
each tube is filled with a sample of homogeneous vol- 
ume fraction, but in practice the volume fraction varies 
slightly from region to region. Additionally, defining 
volume fractions in confined sample chambers is a little 
problematic as the concentration is inherently smaller 



at the walls simply due to packing constraints [57] even 
if the interparticle spacing is spatially homogeneous. 
To define volume fraction, we integrate the r(9) data 
described in the previous paragraph to determine the 
cross sectional area (adding on as to determine the 
physical wall boundary). The length of the observed 
region is known, so therefore we know the volume 
Vtot of the tube that is imaged. Likewise we know the 
numbers of the small and large particles N$ and Nl 
that are in the image, and so the volume fraction can 
be determined from (NgVs + NlVl) /Vtot in terms of 
the individual particle volumes Vs and Vj,. Vs ~ and 
likewise for Vl, so 1 % uncertainties in the particle radii 
lead to 3% uncertainties of the volume fraction. Since 
each particle radius is uncertain, we have an overall 
systematic volume fraction uncertainty of at least 5% 
[49]. There is also some uncertainty between samples 
as the different samples are observed to have different 
number ratios of small and large particles (see Table I), 
and so errors in small and large particle radii will affect 
the different volume fraction calculations in different 
amounts. Encouragingly, visual inspection of the images 
suggests that the calculated volume fractions listed in 
Table I are at least qualitatively in correct order. Samples 
with volume fractions within 0.02 of each other appear 
visually to be the same volume fraction, and samples 
with greater differences are visually distinct. 
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